* 88-1 1169 

CONSTITUTIVE MODELING OF SUPERALLOY SINGLE CRYSTALS 
WITH VERIFICATION TESTING* 


Eric Jordan 

University of Connecticut 


Kevin P. Walker 

Engineering Science Software, Inc. 


INTRODUCTION 

The prediction of component life is an important activity in the development of 
gas turbine engines. Such predictions require accurate component stress analysis. In 
cases involving inelastic strain at elevated temperatures these analyses are 
performed numerically and use constitutive equations as an essential part of the 
analysis. The goal of the research described in this paper is the development of 
constitutive equations to describe the elevated temperature stress— strain behavior of 
single crystal turbine blade alloys. The program involves both the development of a 
suitable model and verification of the model through elevated temperature 
tension-torsion testing. A constitutive model is derived from postulated constitutive 
behavior on individual crystallographic slip systems. The behavior of the the entire 
single crystal is then arrived at by summing up the slip on all the operative 
crystallographic slip systems. This type of formulation has a number of important 
advantages, including the prediction of orientation dependence and the ability to 
directly represent the constitutive behavior in terms which metallurgists use in 
describing the micromechanisms. Along with the model development is a biaxial test 
program that will generate data to determine the model constants, provide insight 
into the qualitative aspects of the material behavior and provide data to verify the 
usefulness of the constitutive model. In the remainder of the report the model will 
be briefly described followed by the experimental set-up and some of the 
experimental findings to date. Finally, several simulations of the experimental data 
using the model will be presented. 

BASIC MODEL FORMULATION 

The inelastic anisotropic behavior of the single crystal superalloy is derived from 
a crystallographic slip model. In nickel— base superalloy single crystals 

crystallographic slip occurs on both the octahedral and cube planes, and it is 
necessary to incorporate both of these slip systems in the model. Indeed, as will be 
discussed in the simulation section, the anisotropic constitutive behavior of single 
crystal superalloys such as PWA 1480 cannot be modelled using only the octahedral 
slip system. 

The constitutive behavior is assumed to be governed by the following set of 
equations: 

*Work done under NASA Grant NAG3-512. 
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where Q . . is the matrix of direction cosines between the global and crystallographic 
axes, and the inelastic strain rate is given by the expression: 
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where 7r r and r r are the resolved shear stress components on the octahedral and 
cube planes, respectively, whilst k t represents the resolved shear stress component in 
the <112> direction on the octahedral plane and produces a yield stress 
asymmetry. The term <j> x is the resolved shear stress on the cube plane in the 

direction of the octahedral slip direction and this term is also involved in the yield 
stress asymmetry. The terms a T and 0 r represent the back stress components on 
the octahedral and cube planes, and K T and L T are the corresponding drag stress 
state variables. The constant, q, in the expression for the drag stress evolution 
equation allows for latent hardening in the single crystal slip systems. At 871C, no 
latent hardening has been detected, and we have accordingly set q= 1. 

The material constants for the simulations were taken to have the following 
values: 

p 1 =4.509xl0 7 , p 2 =645, p s =9.8xl0 -6 , p 4 =0, p 5 =0, K t °= 25775, m=1.49, p=6.69, 5=1, 
p x - 0, 77j=0, g=l, p 6 =1.506xl0 7 , p^lOOO, p^^xll? 6 , I°=71152, n=1.05, t= 7.09, 
/Jg=0, 7? 2 =0, * 2 = 0 , «=0, Z> 1111 =31.37xl0 6 , U U 22=20.84 x10 6 ) £> 1212 =21.26x10 6 . 

The foregoing material constants were determined with the COPES— CONMIN 
general purpose optimization code, which used about one hour of computer time on 
an IBM PC— XT with an 8087 math coprocessor chip. This speed is due to the fact 
that simplified formulae are available for the inelastic strain rate in the <001 > 
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direction for octahedral and cube slip. 


As shown in the simulations it is necessary to include both octahedral and cube 
slip in the theoretical formulation in order to capture the tension— torsion biaxial 
behavior of single crystal superalloys such as PWA 1480. 

EXPERIMENTS AND RESULTS 

The determination of model constants and model verification is being carried 
out using an elevated temperature (87 1C) tension— torsion servohydraulic test 
set-up developed as part of a fatigue research program (NASA NAG3— 160). The 
maximum temperature capability of the system was increased during this program. 
The basic servohydraulic is a 2500 N— M , 250 KN servohydraulic tension— torsion 
machine differing from the commercial machines only by virtue of a University of 
Connecticut developed load frame made from a large die set. Force and torque are 
measured with a conventional strain gage load— torque cell while strain is measured 
by a University of Connecticut designed extensometer that utilizes high temperature 
noncontacting displacement probes. This extensometer, shown in Fig.l, allows 
testing in excess of 871C' without any need to cool the extensometer. The machine 
is computer controlled using a DEC LSI— 1123 microcomputer and software 
developed under this grant. The software allows a wide variety of tests to be run 
including tests involving a sudden change in strain rate, known (cf. reference (2)), 
to be useful in assessing the drag stress state variables used in the anisotropic 
constitutive model. The software also allows strain holds to be inserted into a 
strain history which is particularly useful in assessing the back stress variable 
(reference(3)). The material being tested is a high volume (659t) gamma— prime 
single crystal material, PWA 1480, graciously supplied by Pratt & Whitney. 

The tests run are unique, being the first biaxial tension— torsion tests run on 
this type of single crystal alloy. Biaxial testing has the special advantage that 
orientation effects may be studied by varying the stress state as opposed to growing 
specimens in a wide variety of configurations as is usually done for uniaxial testing. 
Early experiments with this material show that at 871C (1600F) it is possible to 
return the material to the same deformation state by applying a few repeated 
cycles of almost any strain history. As a result a large number of different types of 
tests can be run on a single specimen. Some of the more interesting results will 
quickly be presented. Figure 2 shows strain rate dependence in torsion and it is 
qualitatively similar to tension results. When combined with tension tests these 
torsion results provide a powerful and easy way to determine the relative 
proportions of octahedral slip and cube slip, which is vital to the success of the 
proposed model. This will be elaborated in the next section. Many materials are 
known to cyclically harden to a much greater degree in non— proportional loading, 
reference (3), and it is important to determine if PWA 1480 exhibits this common 
behavior. A test series was run in which a torsional loop was taken followed by 
cycling wdth a torsion to tension strain ratio (X) of 1.5 and with tension and torsion 
90 degrees out— of— phase. This type of non— proportional cycling is known to be 
especially effective, reference (3), in producing extra hardening. Figure 3 shows the 
torsional loops before and after this cycling and the loops are so similar as to be 

nearly indistinguishable indicating that there is no extra hardening. Finally, these 
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alloys are known to exhibit tension-compression yield stress asymmetry which is 
predicted by the proposed model, reference (4). This asymmetry should also 
manifest itself in biaxial tests. Figure 4 shows the asymmetry as a function of the 
lambda ratio (X). It is apparent that at 1600F there is little asymmetry except at 
the smallest lambda ratios. Because the torsion loops are very small at these 
lambda ratios. These findings will be rechecked in an additional experiment. If 

confirmed this data will provide a further test of the ability of the proposed model 
to predict asymmetry. Future testing will include other test temperatures and will 
examine the anisotropy of the back stress and drag stress through the use of two 
rate tests and biaxial strain hold tests. 

SIMULATIONS USING THE MODEL 

The tubular specimens are oriented with the cylindrical axes along the <001 > 
crystallographic direction. When pulled in the axial direction the resolved shear 

stresses on the crystallographic cube planes are zero and only octahedreal slip 
occurs. Material constants for the octahedral slip component of the constitutive 
model were most conveniently derived from the uniaxial data obtained at Pratt Sc 
Whitney. Figure 5 shows the experimental uniaxial data for PWA 1480 at 871 C at 
strain rates of 10 -3 , lO"" 1 and 10 -5 per second in the <001> direction and Fig.6 
shows the corresponding loops predicted with the constitutive model which agree 
well with the data. A check was made to ensure that the Pratt Sc Whitney 

uniaxial data, obtained from a different heat of material and with a bar instead of 
a tubular specimen, was comparable to the data collected in this program. The 
loops generated at Pratt Sc Whitney and at the University of Connecticut at 0.8 
percent strain range and a strain rate of 10 -3 per second were practically 

indistinguishable. Thus, the material constants obtained from the test at Pratt Sc 
Whitney correctly predict the uniaxial data obtained at the University of 
Connecticut. 

Figure 2 shows torsional hysteresis loops carried out at a strain rate of 10 -3 , 
10“ 4 and 5xl0 -5 per second at a strain range of 0.6%. If only octahedral slip is 
included in the constitutive model the predicted loops, shown in Fig 7, are too 
narrow. Inclusion of cube slip in the model widens the predicted torsional loops and 
the predictions with combined octahedral and cube slip are shown in Fig. 8. 
Comparison of the data in Fig. 2 and the prediction from the combined octahedral 
and cube slip theory in Fig. 8 shows very good agreement. It is worth noting that 
<001 > tension data is especially suitable for the determination of the octahedral 
slip constitutive constants, since cube slip is inoperative for this direction under 
uniaxial loading conditions. Torsion tests will involve both cube slip, if operative, 
and octahedral slip and can be used to determine the cube slip constants once the 
octahedral constants have been determined from the <001 > uniaxial data. 
Comparison of the predicted torsional loops, using octahedral slip only, can be used 
to detect the occurrence of cube slip which is usually found only by laborious slip 
trace studies that are very difficult to conduct at temperatures above 593C (1100F). 

SUMMARY 

A slip system based constitutive model has been formulated for PWA 1480 
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which represents explicitly both octahedral slip and cube slip. Tension— torsion tests 
on single crystal material have been run at 871C and some of the torsional results 
heve been sucessfuly simulated using the model. It was found that the prediction of 
torsion data from tension data was only successful if the correct combination of 
cube and octahedral slip was used. The prediction of torsion data is a useful tool in 
determining which slip modes are active. Future work will involve simulations of 
various biaxial test conditions and will also attempt to simulate tests with variable 
strain rates and strain holds. The results of the initial simulations are very 
encouraging and generation of more challenging test cases in the lab and the 
simulation of these more complex cases is expected to be an interesting exercise of 
the constitutive model development. 

REFERENCES 

1. Krempl, E., RPI Report No. RPI CS 78—5, 1978. 

2. Oytana, H.S. and Sidebottom, O.M., "Constitutive Equations Study in Biaxial 
Stress Experiments," J. of Eng. Mat. and Tech., Trans, of ASME Ser. H, Vol. 104, 
Jan. 82, p. 1. 

3. Lamba, H.S. and Sidebottom, O.M., "Cyclic Plasticity for Non— proportional 
Paths: Part I — Cyclic Hardening Experiments," J. of Eng. Mat. and Tech., Trans, 
of ASME Series H, Vol. 100, Jan. 78, p. 104. 

4. Shah, D. and Lin, L., J. Metals, Vol. 32, No. 8, p. 62, 1980. 


281 


TORSIONAL STRESS IN PSI 



Figure 1: Single crystal specimen with extensometer and induction 

heater during biaxial test at 1600°F. 



TORSIONAL STRAIN IN PERCENT 

Figure 2: EXPERIMENTAL TORSIONAL LOOPS FOR P¥A 1490 AT 1600F 
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YIELD ASYMMETRY 


1.0 



0 OA O't. 6.1 0-4 0.5 <X<* 0.7 O.ft l-o (•( \J% m \.V 1-4 t-5 


G 

l.s o 

1-4 * 

Torsional 


\-%- 



O «.i o.J »•> o-H o.5 «•*- o.T o.» 0.4 | .o 1. 1 |.X- \ |,4 l-S 


284 


STRESS IN PSI STRESS IN PSI 


16DD0D 


96000 


32000 


-32000 


-96000 


- 160000 




i 

Ttl DATA 

h 




u 

wk 

Wk 


ppl 

(dl/A 

Wk 


1 

m 

m 



II 






-1.40 


-a 94 


-0.28 


0.28 


0.84 


1.40 


STRAIN IN PERCENT 

Figure 5: PVA 1480 EXPERIMENTAL HYSTERESIS LOOPS AT 1600F 
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Figure 6: PREDICTED LOOPS FOR PVA 1480 AT 1600F (001 DIRECTION) 
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